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Abstract 

In the year 2003 a paper by G. Lapenta demonstrated that there is 
a "new class of soliton-like solutions for the Grad-Shafranov Equations 
(GSE)". The author determined an appropriate pair of transformations 
of the free fields p (fluid field of hydrodynamical pressure) and B z (z- 
component of magnetic induction field) that leads from the Helmholtz 
Equation to the Non-Linear Schrodinger Equation (NLSE) with cubic non- 
linearity. In the following year (2004) , the work of Lapenta was opposed 
by G.N. Throumoulopoulos et al., who criticized his idea of the field 
transformations as a mathematically incoherent choice; contextually the 
authors suggested a new point of view for this one. In his response, in 
the same year, G. Lapenta carried out numerical simulations that showed 
the existence of solitonic structures in a Magnetohydrodynamical (MHD) 
plasma-context. In the present work I want to demonstrate a critical 
condition for a 'complex' poloidal flux function Hf p in a plane framework 
(x, y) that leads to a class of pseudo-general NLSEs which establishes the 
validity of both G. Lapenta and G.N. Throumoulopoulos et al. choices for 
the field transformation functionals. 

*( permanent e-mail address: michele.romeo.mr@gmail.com ) 
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1 Introduction 



In theoretical models in which the accretion dynamics in magne- 
tohydrodynamical plasmas are explored ,i.e. plasmas that obey at 
the various manetohydrodynamical models ,with good approxima- 
tion, in Plasma Physics, like resistive model, Hall-model and others 
mixed and dominant models which include the General-Relativistic 
Theory, we can see that the astrophysical plasmas develop them- 
selves around compact objects, like neutron stars and black holes 
usually; in such theoretical models in which the formations of anti- 
parallel plasma-jets are often examined in the context of their singu- 
lar morphology, there are many references at well defined character- 
istic mathematical structures that could explain the shape of some 
sub-formation inside the accretion and ejection plasma mechanics. 

These mathematical models result very interesting for both their 
internal coherence and physical plausibility, because they can impose 
an indiscutible mathematical link between the fundamental theory 
of the above astrophysical phenomena and the physical charateris- 
tics of the singular plasma formations in the above star structures, 
without affecting the essential set of solutions in the base phisical 
models that explain very well yet the most considerable star dynam- 
ics in the astrophysical plasma context. 

For example, there are many mathematical singular structures of- 
ten investigated because they match at several characteristic cluster 
formations inside the accretion disc and the plasma jets; these clus- 
ter formations are well localized physical structures, with magnetic 
induction field B and electric field E topologies that take place in 
the surrounding space and that donate, with the necessary hydro- 
dynamical force fields, a precise shape at these formations, which 
are therefore confined in limited space areas in the curved space 
around the concerned compact object (this situation exist because 
the gravitational field densities near these astrophysical objects are 
relevant and a full relativistic treatment in any possible theoretical 
model must be considered) and that show in many cases a peculiar 
space periodicity in a well defined global space configuration when 
the magnetohydrodynamical mechanics are in steady state condi- 
tion, or in an equilibrium state. 

These localized plasma formation are indentifiable, for example, 
in the space development of a plasma jet which derives from a com- 
pact star and it stretch itself from the central object to the far space 
areas by a series of quasi-spherical plasma structures named 'plasma 
bubbles' P,[2] (in this case the jet acquired a singular morphology 
known by the name of 'knotty structure') . 

In other cases, similar structures are localized in the accretion 
discs around same compact stars, leading to the formation of related 
structures morphologically, in which the characteristic 'toroidal' sha- 
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pe is linked to the spiral motions in the disc, which response to a 
specifical magnetohydrodynamical equations set of equlibrium state 
that persist on a long time-scale, meaning in this way the same 
steady-state configuration. In this regard, localized structures in 
the dynamics of accretion discs in an astrophysical context have 
been identified analytically and 'experimentally', by several numer- 
ical simulations, from Petviashvili et al. [3],[4j in the 80s as a solution 
to the problem of the Grad-Shafranov Equation (GSE), a funda- 
mental non-linear equilibrium equation in Plasma Physics, by the 
context of a general problem of two-dimensional dynamical plasma 
equilibrium. 

Indeed the Grad-Shafranov Equation \j is a 2nd-order elliptic 
PDE, from the mathematical point of view, which responds to a 
plasma equlibrium problem in a two-dimensional framework {x,y) 
[H] and it's derived by the Euler Equation (EE) in the ideal limit 
(ideal Magnetohydrodynamics) and the Ampere's Law (i.e. 'source 
equation' in the Maxwell's set equations) with neglect the 'displace- 
ment current', because the non-relativistic limit (v << c) or the 
quasi-classic relativistic limit (it takes into account by a 'pseudo- 
newtonian' gravitational potential, like Paczynski-Wiita potential) 
of the gravity potential term in the Euler Equation, which derived 
for the concerned hydrodynamical problem (i.e. accretion disc or 
plasma jet in the astrophysical objects framework). 

This equation establish therefore a 'non-linear problem' for a 
certain scalar 'flux function' 9 which defines the induction field B 
in the ordinary space (x,y, z) around the compact star and which 
have been yet analized from Petviashvili et al. [3] in an analytical work 
that demonstrated the existence of solitary toroidal structures in the 
validity framework of the GSE; 

furthermore, the GSE for this flux function 9 establish the nec- 
essary conditions (i.e. a 'critic conditions') to which the 'parameters 
space' and the free fields in the fluid equations must obey so that 
the equilibrium scenario is respected in the steady-state dynamics. 

At this point, it's indispensable to watch that in many non-linear 
problems of Mathematical-Physics that usually develop some sets of 
expressions containing PDEs, like non linear wave equations, we are 
led to consider, by an appropriate choice of a function transforma- 
tions set and a 'reduction of potentials', an equivalent and more gen- 
eral non-linear differential problem that represents the 'root differ- 
ential problem' for a more general category of non-linear problems. 
From the point of view of this one, the more considerable non-linear 
equation that is present in many problems of Mathematical-Physics 
regarding non-linear wave phenomena that involve PDEs is obvi- 

: ( H. Grad and H. Rubin (1958), V.D. Shafranov (1966) ) 
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OUsly the well known ' Non-Linear Schrodinger Equation' (NLSE), an 

equation that is often present, for example, in non-linear optical 
problems [TT] . 

You know the best known solution scenario of this non-linear 
equation is a ' solitonic solution' or ' 'soliton solution', where a 'soli- 
ton' is the well known special solitary wave solution of several non- 
linear dispersive wave equations; the main characteristics of this 
wave solution are its singular mathematical structure and its phys- 
ical behavior, which synthesize them primarily in their peculiarity 
of localized wave-structure. 

It seems therefore that the GSE, a non-linear (elliptic) equation, 
can be lead to the more general form of the NLSE for the station- 
ary case (steady-state for a plasma in an equlibrium state) by an 
appropriate choice of 'transformation functionals' for the free fields 
inside it; this consideration is furthermore reasoned by the fact that 
the peculiar plasma knots in the accretion discs and plasma jets 
cited above could be justified by an appropriate 're-interpretation' 
of these formations as periodic solitonic structures with justified 
contour conditions. 

At this point, if such appropriate choice of a transformations set 
exist and it is coherent from the mathematical point of view, there 
are no reasons for to refuse such a possibility. Such a choice for 
turning GSE in a NLSE with a specific cubic non-linearity, i.e. the 
NLSE for the Kerr effect in a non-linear optical mean, was made 
recently by G.Lapenta [|], who was demonstrated how it's possible to 
obtain a pair of appropriate transformations in a two-dimensional 
framework (x,y), like that in the equatorial plane of an accretion 
disc, which lead from a general GSE to a cubic NLSE by a two- 
dimensional Helmholtz Equation (2dHE) [7], which is important in 
the non-linear optical problems and it is related to many problems 
for steady-state oscillations (mechanical, acoustical, thermal, elec- 
tromagnetic); but in the following year, by a purely mathematical 
work, Lapenta's idea was opposed by Throumoulopoulos G.N., Hizanidis 
K. and Tasso H. [9\, whom criticized one of his pair of free fields trans- 
formations because its supposed mathematical inconsistency. The 
authors proposed furthermore a new pair of correct transformations 
for the same problem. 

In this work both Lapenta G. and Throumoulopoulos G.N. et al. po- 
sitions will be examined and a new and more general choice will 
be proposed by me, who will show that a 'critical condition' exist 
for the flux function 9, so that this general choice be a coherent 
choice, and both choices are not only mathematically coherent but 
the more general choice cited above include any possible set of rea- 
sonable transformations, leading in this way to a more desiderable 
pseudo-general form of the NLSE [6]. 

All the analytical considerations will be made for the case of 
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knotty structures in the accretion discs as in Lapentas work, therefore 
in the equatorial plane or in the (x,y) framework; in this case z is a 
negligible variable. 

A brief mention to the 'moment method' in general NLSE exam- 
ined in the paper of Garcia-Ripoll, J.J. and Perez-Gartia, V.M. 6J will be 
done in the next paragraph. 

2 From GSEs to a pseudo-general form of the 
NLSE - stationary case 

We derive in the next sections a general choice for the transformation 
functionals of the free fields in the GSEs which leads to a named 
pseudo-general form of the NLSE in the stationary case (steady- 
state case), showing that this choice include both Lapenta G. and 
Throumoulopoulos G.N. et al. choices and putting in evidence that the 
same choices are compatible with one another. 

2.1 Grad-Shafranov Equations set 

The GSEs set, in the more general form for a fully ionized plasma 
(n e = Zrii), can be derive by both the Euler Equations in the ideal 
limit for the ions and the electrons in a two-fluid treatment and the 
Ampere's Law; in this respect, we take the two-fluid equations cited 
in Lighthill, M.J. [13] 



n e m e — v e + v e • Vv e = ~Vp e - M - n e e(E + v e x B) (2.1) 



n.rrii — v, + v; • VvJ = -Vpi + M + rnZe(E + v s ; x B) (2.2) 



where M is the rate of loss of electron momentum, per unit volume, 
by collisions with ions; neglecting electrons and ions inertia (non- 
inertial approximation) and adding both corresponding sides of (2.1) 
and (2.2), we obtain 

= -Vpe - Vpi - n e e(E + v e x B) + n t Ze(E + v 4 x B) = 





Vp + n e e(-Vi — v e ) x B 



Vp + JxB 



which lead to 



Vp = J x B 



(2.3) 



being the current density 



J = n e e(v s: - v e ) 
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At this point the calculus is the same of Lapentas work. Assuming a 
plane framework (x, y) because the analysis is made on the equatorial 

plane of an accretion disc,we can take as a particular solution of (2.3) 

m 

B=ixVf + B 2 2 (2.4) 
The Ampere's Law lead by the (2.4) to 

J = V x B = z(V • V*) - W(V-z) + (V* • V)z + 
- (z ■ V) V* + B z (V x z) - z x VB Z 

which delivers, being z a negligible variable, 

J = zV 2 * — z x VB Z (2.5) 

If we use the current density expression (2.5) and the solution (2.4) 
in (2.3), we obtain the equation 

Vp = (zV 2 * - z x VB Z ) x (z x W + B z z) = 

= zV 2 * x (z x V*) - (z X VB Z ) x (z x V*) - (z x VB Z ) x B z z = 
= V 2 * (zz • W - • z) - [(z x VB Z ■ V*) z - (z x VB Z ■ z) W] + 
+B Z (zz • VB Z - VB z z ■ z) 

which delivers 

Vp + V 2 + {zx VB Z ■ W) z + B Z VB Z = (2.6) 

Observing that equation (2.6) is made of two parts linearly indepen- 
dent among them, from this one it follows the equations system 

(z x VB Z ■ V*)z = 

Vp + V 2 + B Z VB Z = 
which represents the 'equivalent GSEs set', or 

WxVB z -z = (2.7) 

Vp + V 2 yV$> + B Z VB Z = (2.8) 
where the equation (2.7) derives from 

z x VB Z • V* = (2.9) 

The equations (2.7) and (2.8) are the vectorial-scalar form system 
of the GSE for the ideal magnetohydrodynamical equilibria in a 
two-dimensional framework and they lead to a fundamental 2nd- 
order PDEs system for searching the steady-state solution of the 
stationary problem for a fully ionized plasma (astrophysical plasma). 
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2.1.1 Canonical form of the GSE 

We can observe that the canonical form of the GSE derives from 
the system (2. 7), (2. 8) as it follows: by the equation (2.7) we have 
necessarily 

Vf xVB z =0 V{x,y) E Dom(^)nDom(B z ) (2.10) 

because V* x VB Z is on the z axis; this one delivers a linear depen- 
dence between V* and VB Z as it follows 

Vf = a(*,5,)VB z (2.11) 

where we suppose that a is a rational real functional of * and B z 
defined in Dom(*) n Dom(B z ); therefore we obtain by the (2.8) 

Vp + V 2 *a(*, B Z )\7B Z + B Z VB Z = 

which delivers 

Vp+(V 2 ^a{^>,B z ) + B z ) VB Z =0 (2.12) 

this one means that VB Z and Vp are linearly dipendent among them, 
or they are co-axial vectors; furthermore, by applying V*x to this 
one, we obtain 

V* x Vp = V(ar, y) e (£>om(*) n Dom(p)) (2.13) 

which means that V* and Vp are necessarily linearly dependent as 
in (2.10); thus determines 

V* = /?(*, p)Vp (2.14) 

where we suppose that is a rational real functional of * and p 
defined in Dom(^) n Dom(p); now, we can observe that (2.14) lead to 
the differential equation 

= V* - ,p)Vp = 

which is 

= D(V,p) V(x,y) e (Dom(y ) n Dom(p)) (2.15) 

where we suppose that D is a real differential functional dependent 
on (3 or 

D(^,p)(x,y) e R V(x,y) G flom($) n Dom(p) 
a general integral of (2.15) can be written as 

p = Fp(V) V(.x, y) eftp C (Dom(V) n Dom(pj) (2.16) 

where is a real functional dependent on fi what is defined in tyg; 
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If we take the equation (2.11) instead, we have a differential equa- 
tion equal to 

= V* - a(tf , B Z )VB Z = B z ) 

which is 

= £>(tf,.B z ) V(i,i/)e(Dom(*)nDom(B z )) (2.17) 

where we suppose, as above, that D is a differential functional de- 
pendent on a; a general integral of (2.17) therefore can be written 

as 

B z = V(x, y) e Q C (Dom($) n Dom{B z )) (2.18) 

where F Q is a real functional dependent on a and it's defined in Q a ; 

now, if S ^ is the unit vector for both VB Z and Vp, from (2.12) 
it follows 

||Vp||6+(V 2 *a(*,B z ) + B z ) ||VB,||6 = 
[\\\7p\\ + (\7Ha(y,B z )+B z )\\\7B z \\]b = 
\\Vp\\ + (vHa(y,B z )+B z ) \\VB Z \\ = 

which lead to 

V 2 * = (- ||Vp|| HV^ir 1 - B z ) a(*, -Bz) -1 (2.19) 

and taking in account (2.16) and (2.18) in this one, we obtain finally 
the GSE in canonical form 

V 2 vf= (-||VJJ 3 (*)||||VJ'„(*)|r 1 -J , „(*))a(*,F a (*))- 1 (2.20) 

in which we can recognize an Helmholtz Equation. In this one we 
suppose that 

, \ im , , ^TfM\ eR ' vfav)en a nn p :\\VF a (*)\\ = o 

(x,y)-*(x,y) \\\/l< a [W)\\ 

because the flux function * has not any pole inside the equilibria 
domain Q a n Q/3. 

2.2 A general choice in the complex plane for the free 
fields 

A pseudo-general form of the NLSE show at first view a more great 
complexity than an ordinary cubic NLSE cited in Lapenta 's paper 
but on the other hand its mathematical morphology is affected by a 
more general set of solutions, which can represent a wide-spectrum 
of possibilities for the research of explanations in the physical frame- 
work for the above equilibrium structures. Furthermore, the general 
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form of the NLSE cited below includes both Lapenta choice [8] and 
Throumoulopoulos et al. choice [9] (this last leads to a pseudo-cubic 
NLSE); in this way, several specialized choices set for a same phys- 
ical problem can be put in a more general mathematical solving 
context from the point of view of the related non-linear problem. In 
this case, such a strategy therefore permits us to define a valid vay 
for to find a more general possible solving method for a GSEs set in 
object, working in the framework of the solitonic solutions as in the 
interesting idea of Lapenta. 

A pseudo-general form of the NLSE, in the time-dependent case, 
can be represented as it follows [6] 

A 3 ^ + g (\V\ 2 ,t) * + ia(|*| 2 ,t) * (2.21) 

where * is a complex function, A 3 is the Laplace operator in the 
(x,y,z) space framework and both g and a are supposed to be time- 
dependent potential complex functionals of |*| 2 . It's remarkable the 
fact that this equation is named 'pseudo-general' because the func- 
tions g and a are complex- valued functionals instead that real- valued 
functionals as is usually. However, it's clear that (2.21) can be put 
'always' in a form in which g and a lead to a pair of real-valued 
functionals. Indeed 

.9* + icr* = (g R + i gi ) ty + i{a R + ia t ) * = 

where g and a are obviously real functionals. In the stationary case, 
this equation leads to the steady-state form 

A 3 *+ (.9 (|*| 2 ) + ia (l*| 2 )) * = ( 2 -22) 

which is, as we see, a more general form of the Helmholtz Equation; 
it can be viewed as [7] 

(A 3 W(|*| 2 ))vK = o (2.23) 

where 

^ 2 (|*| 2 ) =g(\M 2 )+^(\M 2 ) (2-24) 

is a complex-valued function. This non-linear PDE therefore repre- 
sents the fundamental link between the GSE and the pseudo-general 
NLSE in the stationary case for a three-dimensional equilibrium 
problem. In the two-dimensional framework (x,y), that is the frame- 
work in which we analize the localized accretion plasma structures, 
we have therefore the differential equation 

(A 2 +ro 2 (|*| 2 )) * = that is V 2 * = -w 2 (|*| 2 ) * (2.25) 
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At this point, we observe that an appropriate choice for the free 
fields in the GSEs can permit the translating of such equilibrium 
equations in a certain form of the NLSE, as the general stationary 
form (2.22), which could be an excellent solution for re-interpretate 
localized structures in the accretion discs as 'solitonic structures', 
going as viewed by the Helmholtz Equation. Indeed we'll see that 
a general transformations set for the free fields p and B z can lead 
from the equations (2.7) and (2.8) to the equation (2.25) by an 
appropriate expression for the w 2 function on the same p and B z . It's 
important to observe that the 'analiticity' of the complex function 
* is not required for solving the Helmholtz Equation in general [TJ. 

Furthermore, I cite that the 'moment method' [6| can be used as 
a valid approximation analytical way for to solve a wide family of 
non-linear wave equations of NLSE type in its pseudo-general form; 
this method have been developed for n-dimensional cases in general 
by Gartia-Ripoll, J.J. and Perez-Garcia, V.M.. So, it's remarkable that 
by this mathematical approximation strategy it could be possible 
to identify a good set of physically coherent solutions of the GSE 
for the magnetohydrodynamical equilibrium problem, i.e. solitonic 
solutions as we would like. 

2.2.1 Connection between GSEs and a pseudo-general form of the 
NLSE 

Now, we consider the expressions (2.16) and (2.18) in the intersection 
of the B z and p domains 

P = F^) B z = F a (^) V(x,y) eflanflp 

and replacing the free fields p and B z in the equation (2.8) with this 
relations we obtain 

Vf>(*) + V 2 *V* + F Q (*)VF a (*) =0 

obviously for the linear dependence of Vp + B Z VB Z by V 2 ^fV^f this 
one delivers 

Vi^(*) +F Q (*)V J F a (*) = fc(tf) V* = -V 2 WW (2.26) 

or 

(V 2 * + k (*)) V* = (2.27) 
which is V 2 * + fc(*) = V(x, y) e fl a n ftp 

that is exactly the Helmholtz Equation, in which we suppose that 
k is a generic functional of if we compare the second equation 
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in (2.27) with the Helmholtz Equation (2.25) and we consider the 
expression (2.24), then we obtain necessarily 

fc (*) = ro 2 (|lf) = (.g^*! 2 ) +ia(|*| 2 )) * V^t/je^n^ (2.28) 

which indicates that k is a complex functional and * is clearly 
a complex function by the equation (2.27), because by this one it 
derives 

V 2 * = -fc(tf) = - [g (|*| 2 ) +ia (|*| 2 )) * 
or *- 1 V 2 * = - (.g (J*| 2 ) +ia (|*| 2 )) with g, a ^ (2.29) 

At this point, taking into account the Lapenta and Throumoulopoulos 
et aVs works [S], [9] and starting from their fields choices, I found a 
well posed pair of transformation functionals like this 

B Z VB Z = ivB 2 = £)aylV +1 V¥ (2.30) 

i,j 

Vp = ^6yF^ +1 V* (2.31) 

where * is the complex conjugate function related to the complex 
flux function ^ and a iJ: are complex constants; we observe that 

(2.30) and (2.31) are respectively a term proportional to the mag- 
netic force density and the pure hydrodynamical pressure force den- 
sity, which are two 'real' physical variables; this is for me a funda- 
mental physical condition and it's important that it doesn't affect 
the complexity of the flux function ^; we will explore this 'reality 
condition' in the next section as a 'validity condition' for (2.30) and 

(2.31) . From these transformations and the first equation in (2.26) 
it derives 

Vp + B Z VB Z = (aij + 6; 3 ) * = k (*) W (2.32) 

hi 

which together with (2.28) delivers 

fc W = ^K+M* l ^ +1 = (2-33) 

= (.g(l*l 2 )+^(l*l 2 ))* V(ar,y) eO a nfl^ 

Now, it's interesting to observe that the fields transformations (2.30), 
(2.31) lead from equations (2. 7), (2. 8) to the pseudo-general form 
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(2.25) of the NLSE in a space framework (x,y) by the equation 

Et^***' = g (\*\ 2 ) +t* (\*\ 2 ) V(x,y)en a nfi (2.34) 

where 7^ = + b io e C and Im(*)cC 

Because of the second equation in (2.29), we will proof below that, 
if * is a complex function, X/y^-'Sr*- 7 ' can be put in the form 

. 9 (|*| 2 )+«7(J*| 2 ) 

where both g and a are supposed to be complex functionals as above. 

In this regard, we consider the expression at left hand side of 
(2.34) in Q a n ftp domain 

E^,* V (2-35) 
and if we develop the sum, we obtain the next identities chain 
Y.':^' xl ' J = E" ' • >l ' ' + E^ - **' + E^* 2 *' + - = 

i,3 3 3 3 

Joo+7iol v I / | 2 + 72o(l*| 2 ) 
+ (7oi+7nl*| 2 +72i (l* 

= E^o (W) *° + E^i 0i 2i ) + E^2 (W) * 

i i 

= E(E^(i*i 2i ))^ = E/, (i*i 2 ) 

which deliver 

E7«*** j ' = E^(i*i 2 )* j ' (2 - 36) 

i,j 3 

where 

.^(i*i 2 ) =E^-(i*r) ( 2 -37) 

i 

is a complex polynomial functional on |*| 21 . Now, we 'impose' the 
complexity of the flux function * (i m (*) c C) and we write therefore 

*(x,y)=u(x,y)+iv(x,y) V(s, y) € £l a n fi/j (2.38) 



,2> 2 



* + ... = 



2 + ...= 
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where u and v are generic real- valued functions in the two-dimensional 
framework (x,y); from this one and (2.36) we obtain 

x;7«* < * j '=E^(i*i 2 ) (u+iv)j = (2 - 39) 

i,j 3 



k=0 x ' j k=0 



in which we used the 'Newton binomial formula' and ( Jk = ( 3 k ) are 
the binomial coefficients. 

At this point, we fix the maximum level N of the polynomial sum 
in (2.39) and we observe that 



AT j N j 



E^(i*i 2 )E^ ui_fc ( iv ) fc = ED 0*i 2 ) ^u j - k (iv) k = 

3=0 k=0 j=0 k=0 

= EI>* (1*1 V~ fc(it,)fc ( 2 - 4 °) 

where „ jfc (|*| 2 ) = f 3 (|*| 2 ) ( Jk = (f)E^ (j^) ( 2 - 41 ) 

is a new complex polynomial functional; from (2.40) therefore it 
derives 

N j 

EE^(l*| 2 ) ^- k {iv) k = 

j=0 k=0 

N N-l N-2 

= E%> (i*i 2 ) u l (ivf+j2v { i +m (i*i 2 ) ^m 1 +E^+2) 2 (i*i 2 ) «'m 2 +... 

;=o z=o z=o 

that is 



A72 /N-2h \ N/2-1 (N-{2h+\) 



E E m m (l*l 2 ) «' (-) 2h + E E ^(a/H-i) (l*l 2 ) «' (<«) 

ft=0 V Z=0 / fc=0 \ i=0 / 

if iV is 'even'; in clear complex form this one deliver, 

E il " = ( 2 - 42 ) 



2/i+l 



E E^)(i*i a ) tt ' (- 1 ) h « afc + 

'i=0 \ 1=0 ) 

f-1 (N-(2h+\) \ 

*E E (!*!>' (-!)^ : 



+ I 



2h+l 
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and if it's compared with equation (2.34), taking into account (2.41), 
it delivers finally the next expressions for the g and a functionals 

V(ar, y) 6 0„n0^ 

f /N-2h \ 

. 9 (l*| 2 ,^) = £ I £ (l*l a ) (-i)V> = 

h=0 \ (=0 / 
f iV-2h / , N 

=EEEuV«K)«'(- 1 )*« tt ( 2 - 43a ) 

2 h=o ;=o v y 

f-1 /N-(2h+l) \ 

,(|*| 2 ,U,,)=£ £ ^ + l)(l*l>' (-l)^ 2ft+1 = 

/i=0 \ 2=0 / 

#iV-(2ft + l) 

= EE E L+J^O*! 24 )"'^ 1 ) ^ 1 ( 2 - 44a ) 

which are clearly two complex functionals on the square function 
l^l 2 = and on the u and v functions; they are complex because 
the presence of the complex coefficients -f u in the above expressions. 

If the level of the polynomial sum N is odd, we have similarly the 
complex functionals 

N — 1 

2 (N-2h \ 

g(\*\\u,v)= W E^) (\*\ 2 )u l )(-l) h v 2h = 

h=0 \ 2=0 I 

N — 1 

—2~N-2h , . \ 

= E E E ( 2 J7«(l*| W )t.'(-l) h « afc (2.43b) 
i h=0 i=o ^ ' 

(N-(2h+l) \ 

ff (i*i 2 , = £ e »&(*h-i) (i*i a ) «' (-i)^ 2ft+i = 

/i=o y 2=0 J 

^l N -(2h+l) 

= EE E L+iJ^O*! 9 *)^^^ 1 ( 2 - 44b ) 

i h=0 1=0 ^ ' 

We observe that in these expressions the maximum level of poly- 
nomial sum on i-index is not fixed, leaving in this way the choice 
of dependence of the above functionals on the even powers of the 
function |*| 2 completely free. Such a dependence can be imposed on 
the base of mathematical or physical criteria which are well defined 
if they are related to a specific physical problem. 
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Notation 1 it's clear that the general choice (2. 30), (2. 31) includes 
both Lapenta JEj and Throumoulopoulos et al. JM/ choices; indeed for the 
first choice we have 

a m =al, a.y=0 Vi,j > (2.45a) 

b n =a 2 , bij = V(i,j)^(l,l) (2.45b) 
while for the second choice we have 

aoo = &o> a ij =0 Vi, j > (2.46a) 

&03 = «o, &tf=Q V(*,j)^(0,3) (2.46b) 

3 Validity conditions in the mathematical and 
physical frameworks 

In this section we explore the validity conditions for the general 
choice (2. 30), (2. 31) from the physical and mathematical points of 
view. We observe that these conditions are named 'critical condi- 
tions' 1 in the present paper because they are 'necessary' conditions 
for the acceptability, or internal coherence, of the proposed mathe- 
matical models, which are analized in respect of a related reasonable 
link between the mathematical framework and the physical one. 
Furthermore, the analysis of these theoretical positions will take 
into account necessary and sufficient conditions too for the valid- 
ity of unique solutions for the magnetohydrodynamical equilibrium 
problem proposed. 

3.1 On a critical condition for the general choice: mathe- 
matical aspects 

We consider the general choice for the free fields p and B z in the 
GSEs set (2.7), (2.8) 

B Z VB Z = \*JB 2 Z = ^.J^'-'Vt (3.1) 

i,3 

Vp = 2&ylV +1 V¥ (3.2) 
with , b i:j e C 

and for to establish a validity condition for the internal mathemati- 
cal coherence of these transformation functionals, we apply the op- 
erator Vx at the both sides of (3.1), (3. 2), as in the Throumoulopoulos 
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et al. work [9j, taking into account that they are 'similar' except for 
the coefficients sets {aij} ij( . N and {hj} ij& N ; for this reason, it will be 
sufficient therefore to consider for both relations the term 

5^cy*V +1 W (3.3) 

hj 

with Qj = ciij or bij as appropriate, and to reduce in normal form the 
expression 

Vx^Cj/j&V +1 Vtf (3.4) 

hj 

for to identify a validity condition which is unique for both field 
transformations. At this point the calculus is quite simple: starting 
from the expression (3.4) we obtain for (3.1) and (3.2) 

Vx^cy*V +1 V* = V(x,y) en a nn fJ (3.5) 

since 

VxVp = Vx Q VB *) = 

for the 'scalar' functions p and B\. If we proceed in the calculus we 
obtain 

= Vx^c tJ f$ J+1 V* = ^c i3 Vx (¥V +1 V*) 

i,3 hi 

for the linearity of vectorial product; taking in account the vectorial 
identity 

Vx (a~a) — aVx~a — ~axS7a with a as a scalar variable, 
we have therefore (we remember that * is a scalar function) 

i,3 hj 

= - ^Cy W x V = - ^c I3 V* x ^' +1 vf +f'v* J+1 

hi i,j 

= - ^ ClJ V* x (V+W^V* + 1)^V* 
= - ^c^'+W^W xVf = -J |*| 2 ) V* x V* 



or 



J(f,|$| 2 V$xV* = V(x, y) e n a n n (3.6) 
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where 

j(*,|*| 2 )= X'V"' ' ;7ir ; ( 3 - 7 ) 

is a complex polynomial functional of finite degree on * and |*| 2 . As 
we can see, the equation (3.6) has the solutions 

|*| 2 )=0 (3.8a) 

V$xVf = (3.8b) 

The (3.8a) solution functional equation has no solutions for an 'ar- 
bitrary' flux function because the terms inside it are linearly 
independent, or 

Vi^T,Vj^j Jw e C \ {0} : tfV = ujW^ 
indeed, for this reason it derives, taking in account (3.7), 

V*, ${ Cij } iijeN , ca GC: j(*,|*| 2 ) =0 

this means that the functional equation (3.8a) cannot represent a 
necessary condition so the equation (3.5) is verified; at this point 
it's clear that the unique 'critical condition' for the mathematical 
coherence of the general fields choice for an arbitrary flux function 
* is 

V*xVf=0 V(x, y) e tt a n tt fj 

Now, if we impose the complexity of this function (Im (*) c C) and 
we write therefore (see position (2.38)) 

*(ar, y) = u(x, y) + iv(x, y) V(x, y) e Q n dp 

we obtain, starting from the equation (3.8b), 

C = V* x V* = V (u + iv) x V (u - iv) = (Vu + iVv) x (Vu - iVv) = 
= -iVu xVul iVv xVu = -2iVu x Vv 

(where C = 0+io); in this way, the condition (3.8b) is equivalent to 
the condition Vu x Vv = 0, or 

V*xV$ = c o ViixV« = V(x, y) E fl a n Q,p 

We have therefore the 'critical condition' (which is a 'necessary con- 
dition') for both real and imaginary parts of * 

Viixv« = o v(x, y) e n a n n f3 (3.9) 
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so the equation (3.6) is verified; an expansion in a two-dimensional 
framework (x, y) of the V operator finally leads to the critical PDE 
for the real functions u and v 

„ „ ( du^ du^\ ( dv ^ dv^ 
= VuxV«= — x+— y x — x+— y 
\ox ay ) \dx ay 

dudv ^ ^ du dv ^ ^ / du dv du dv\ 

dx dy X y dy dx y X \ dx dy dy dx ) X y 

or the equation 

( du dv du dv 1 

k%-%s}<-»»=° (3.10) 

At this point, it's clear that the differential equation (3.10) is the 
unique necessary condition for an arbitrary complex flux function * 
so the general free fields transformations (3.1), (3. 2) are a mathemat- 
ical coherent choice for the equilibrium problem in object; further- 
more, it's remarkable that this condition doesn't imply the analitic- 
ity of the function ^; indeed, if the analiticity is requested in fl a nftp 
for ty, the fun ctions u and v must be C - differentiable in the domain 
A = (fl a n flfj) \ d(fl a n ftp), or the Cauchy-Riemann conditions must 
be valid in A; these deliver for u 

I / du \ 2 f du 



while for v we have similarly 




(x,y)=0 V(x,y)eA (3.12) 



and it's clear that the differential equations (3.11) and (3.12) are 
specific cases which are included in the differential equation (3.10). 

Usefully, we note that if the complex flux function * is in the 
Gauss form 

^(x,y)=p(x,y)e^ (x ^ V(a;, y) € fl a H Qp (3.13) 

the critical condition (3.9) must return 

VpxV^^O V(x,y)GQ a nnp (3.14) 

Let us observe that if we impose the analiticity of the flux function 
in the equilibria domain (Q a n Qp) , the Helmholtz problem (2.27) 

V 2 * + fc(*) = V(x, y) e n a n Qp 
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deliver the functional equations system 

V 2 * = (3.14.1) 
fe(tf) = V(x,y) e O a DQ (3.14.2) 

these equations destroy obviously the related NLSE-problem but 
they represents also a not-banal question from the closely functional 
point of view, because the equation (3.14.2) must be satisfied only 
in a 'sub-domain' Q a n^ of the entire domain Dom(*); indeed these 
domains are subject to the condition 

(o Q no^) c Dom(-$) 

leaving free in this way the choice of the not-bound portion of ^ 
inside Dom(^) \ (O a nOa). This means that a certain family of \E» 
functions which satisfy the equation (3.14.2) can be identified as 

$ = h (Q a n 0/3) = * V(a;, y) e (O a n 0/3) n Dom (*) 

and it means also that this family must depend only on the pos- 
sible forms of the flux functions inside the complementary domain 
Dom(f) \ (12 Q nfi,3). Now, because the (3.14.1), we can note that 
if (0 Q nO/j) ^ Dom(^), it can be possible to define a family of har- 
monic functions $ inside O^nO^ which satisfy the functional equation 

(3.14.2) ; this functions family is related obviously to the solutions 
set of the equation 

ff (|*| 2 )+«r(V| 2 ) =0 V(x,y) eO Q HO/3 (3.14.3) 

At this point, it's remarkable that this equation necessarily deliv- 
ers the above not-trivial solution for <5, because there is the close 
condition 

(O a n O/3) C .Dom (*) (3.14.4) 

which can be justified by observing that the equilibria domain does 
not concern in general the entire framework of the equatorial plane 
in the accretion disc unlike 

Let us note that in our specific case the solutions set of the 

(3.14.3) is not trivial because the general equation scenario (3.14.3) 
assume the particular form 

g (|*| 2 , u, v) +ia (V| 2 ,«,«)= V(z, y) £ O a n n f3 

related to the expressions (2.43a)-(2.44b). It's clear that this case 
is generally independent of the condition (3.14.4). 

It's important too, in conclusion, to observe that the 'knotty 
structures' can be present only inside a sub-domain O a nO/3 for which 
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the condition (3.14.4) is valid, giving sense in this way to the plasma 
equilibria in the related space regions; hence, it's necessary that the 
contour conditions for the plasma equilibrium-dependent differential 
Helmholtz problem 

V 2 * + (g (V| 2 ) +«j (l*| 2 )) * = 

must be specified on the contours of an appropriate domain O, taking 
for this f2 c Q a n Qp. 

Notation 2 The critical condition (3.9) is identically satisfied for the 
Throumoulopoulos et al. transformations JMj, while for the Lapenta trans- 
formations f8jj we have that the above condition remains valid in its 
general form, leading to a free choice for the real flux functions u 
and v, which in this case are related only to the bond equation (3.10). 
This is because in the second choice is present a term equal to 
while this one there is not in the first choice. Such situation means 
that both choices represent coherent transformation models but the 
Throumoulopoulos et al. choice is mathematically more strong than 
Lapenta choice; from this consideration clearly doesn't derive that the 
choice (2. 45a), (2.45b) is less valid than the choice (2. 46a), (2. 46b). 
It's important instead that this observation reconciles in this way the 
results of the above authors. 

Notation 3 One of the solutions for the equation (3.6) is 

j(*,|*| 2 )=o \t(x,y) e n a nQ 

Let us note that this functional equation has not a trivial solution 
if we consider an appropriate family of flux functions $ as we did 
above in a similar problem; in this way, we could obtain another 
valid critical condition on the ^ function for the general fields trans- 
formations. An idea is to study the character of the functional power 
series 

00 

for an appropriate choice of the constants set {c^-}. . £ N , cy e C, tak- 
ing into account that it must be 

Joo (*,l*l 2 ) =0 v(x,y)en a nn fj 

3.2 On a reality condition for the free fields: physical as- 
pects 

In previous section we talked about a plausible physical validity 
condition for the general choice (3.1), (3. 2), saying that the 'reality' 
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of the terms proportional to the magnetic force density and the pure 
hydrodynamical pressure force density inside it can be a fundamen- 
tal condition from the physical point of view and it's remarkable that 
such condition doesn't affect the complexity of the flux function 
indeed these terms are the two 'real' vectorial physical variables, by 
the free fields transformations, 



VB 2 Z = 2^ay*V +1 W 



Vp = ^6 ij '*V +1 V# 



which must be therefore 

VB 2 z} \7p^R 3 c , where 



{x c : x c = x + iy,Vx,y G 



(3.15) 
(3.16) 

(3.17) 



Now, for to determine this reality condition for the above force den- 
sities, we impose 

Im (VB 2 ) = Im (Vp) = (3.18) 
and taking the polynomial functional (2.34) 

= 9 (l*l 2 ) +i* (l*l 2 ) , 7« e c, v(x, y )en a nn 



for (3.15) and (3.16) we obtain 



VB 



5 (|*| 2 )+ia(|*| 2 ) 
\7p = [g (|*| 2 ) +ia (|*| 2 ) 



(3.19) 



(3.20) 



which are different only for a factor 2. At this point, for the valid- 
ity of both (3.18) relations, it's sufficient obviously to analyze the 
condition rr , , 

Im| 5 M*| 2 J +i<7 M*| 2 J #V*J=0 (3.21) 

For this one, taking in account that 

#(ar, y) = u(x, y) + w(x, y) V(x, y) e Q a n 0/3 

let us calculate in clear complex form the expression 

3 (|tf| 2 ) -M<r (|tf| 2 )] tfV* 

We have therefore 

(g + ia) = (g + ia) (u + iv) V (u + iv) = 
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= [(gu — av) +i (gv + au)] (Vu + «Vu) = 

= [(gu — av) Vu — (gv + au) Vv] + i [(gu — av) Vv + (gv + gu) Vu] 

which delivers for the condition (3.21) 

(gu — av) Vv + (gv + au) Vu = 

or 

Vu = [ aV - 9U ] V(i,|/)e!! a n[!« (3.22) 
(gv + au) 

where the form of the functions a and g is derived by the relations 
(2.43a)-(2.44b); it's clear that (3.22) represents a necessary and suf- 
ficient condition for (3.21) and furthermore it leads to the 'new' 
critical condition 

VuxV» = V(x, y) £ n a n flf) (3.23) 

It's important to observe that this one is the same of the critical 
condition from the mathematical point of view, which can be found 
in (3.9); this means that the critical condition (3.23) ia a 'common 
necessary condition 7 for both the mathematical and physical aspects 
of the coherence for the transformation functionals (3. 15), (3. 16). 

Now, if we note that the generic solution choosed in a plane 
framework (x,y) for the magnetic induction field is (by the (2.4)) 
V(a;, y) £ O a n 

B=zxV* + B z z 

is reasonable to think that it must be for a complex as in the 
conditions (3.18), 

Im(B)=0 V(ar, y) £ Q n ftp (3.24) 
which delivers necessarily 

= Im(z x V* + B z z) = Im(z x V*) = Im (z x (Vu + i\7v)) = 
= Im (z x Vu + izx Vu) = z x Vu 

or 

z x Vu = (3.25) 

this condition leads to below vectorial differential equation by the 
expansion of the operator V 

^ ( dv ^ dv^\ dv ^ dv ^ 

z x ( — x + — y ) = — y - —. x = 

which means 



dx dy ) dx dy 



dv dv 

— = — = or Vu = V(i,rfe!! a n!lfl (3.26) 
ox ay 
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At this point it's indispensable to observe that if the imaginary 
component of the flux function * has a zero gradient in J] a nO^ the 
critical common condition (3.23) is identically verified; so, the only 
position which can respect a reality condition for all fields in the 
equilibrium problem, as viewed, is therefore a unique flux function 
condition for which it must be 

Im(^(x,y)) = v(x,y) = cost V(x, y) £ Ct a n £1/3 (3.27) 

which delivers finally 

ty(x,y) =u(x,y) +iA, A e E, V(a;, y) £ O a n fi/3 (3.28) 



3.3 On a global critical condition for the GS equilibrium 
problem: critical equations system 

Here we briefly explore a further validity condition for the accept- 
ability of the general choice (3.1), (3. 2). This condition regards the 
equation (2.8), which leads to a final 'validity equation' if it's re- 
lated to the common critical condition (3.9) or to the flux function 
condition (3.27), taking into account the relations (2.11) and (2.14). 

If we start from the equation (2.8), we obtain by applying the 
operator Vx to both sides of this one 

C = Vx (Vp + \7 2 ^V^ + B Z VB Z ) = VxVp + VxV 2 *V4' + VxB z \7B z = 
= V 2 *VxV* - V* x VV 2 * + B z VxVB z - VB Z x VB Z = 
= -V* x VV 2 * = VV 2 * x V* 

or 

VV 2 * xV* = O c (3.29) 

This further condition is obviously a necessary condition for the 
Grad-Shafranov equilibrium problem and it can be translate into 
below differential equations set: if * is the complex function *(ir, y) = 
u(x, y) + iv(x,y) V(x, y) £ Q a n flp, we have 

C = VV 2 *xV* = VV 2 (u + iv) xV (u + iv) = 
= (VV 2 u + iVV 2 t)) x (Vu + iVv) = 

= VV 2 u x Vu - VV 2 !) xVd + i (VV 2 u xVul VV 2 u x Vu) 

which delivers the 'critical equations system' 

VV 2 u x Vu - VV 2 u xVd = (3.30a) 
VV 2 u x Vu + VV 2 !) x Vu = (3.30b) 
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Now, if we consider the flux function condition (3.27), we have that 
the second equation is identically verified, while for the first equation 
we obtain 

vv 2 u xVu = o v(x, y) e o Q n n (3.31) 

or 

= VV 2 «xV!i= -x|-y V 2 ux -x+-y = 
\ox dy J \ox dy J 

d „o du _ _ d o du _ _ 

/ d ^ 2 du d 2 du\ ^ 
\dx dy dy dx J 

which delivers the unique differential equation for all critical condi- 
tions 

If we consider instead the common critical condition (3.9) only, 

ViixV« = V(a;, y) G Q n 0,3 
we obtain by the system (3. 30a), (3. 30b) 

VV 2 m x Vu = V(x, y) G Q n 0,9 (3.33) 
which delivers, as in (3.32), the differential equation 

taking into account that the condition (3.9) means 

Vu = / (u, v) \7v 

where / is a real functional on u and v, V{x,y) e O a n 0/3. 

At this point, it's clear that the conditioning equations (3.32) 
and (3.34) represent two global critical conditions for the magneto- 
hydrodynamical equilibrium problem in the pseudo-general NLSE 
framework, for which obviously it's true that 

f d ^ du d ^ 2 du\ q ^ 
\ dx dy dy dx J 

for the bond equations system (3. 30a), (3. 30b). 
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4 Comments 



In conclusion: we have seen in previous sections that for to ob- 
tain a pseudo-general form of the NLSE in the stationary case (see 
(2.22) in the plane framework {x,y)) by the GSEs set (see the sys- 
tem (2. 7), (2. 8)) for a general plasma equilibrium problem in the 
equatorial plane of an accretion disc, or 

V 2 * + (g (V| 2 ) +icr (l*| 2 )) * = 

a set of transformation relations (i.e. the general choice) of the free 
fields in the GSE is imposable in a general form and this is 

B Z VB Z = ivB 2 = £ay*V +1 V* 

Vp = ^2bijW^ j+1 V* with a t3 , bij G C 

where \t is the flux function in the solution (2.4) of the general 
equilibrium problem (2.3) and it's a complex function on the real 
variables x and y 

y) = u(x, y) + iv(x, y) V(x, y) G tt a n Sip 

which is a solution for the Helmholtz problem (2.25) and it's defined 
in a specific local not-banal domain (in the accretion disc frame- 
work) as flanflp] it's clear that the topology of this domain is a very 
fundamental characteristic of the equilibrium problem, because the 
critical condition equation (3.10) must be verified for the mathe- 
matical coherence of the above fields transformations inside it, or 

( dudv du dv 1 , , „ w , , ^ ^ 
{d- X d-y-d-yd- X \ {x > y)= ° VOr.deO.nn, 

which, in the vectorial form, is 

VmxVd = V(.t, y) efl a r\ flp 

and therefore it's remarkable that the analiticity of the flux func- 
tion is not necessary for the validity of the general fields choice. 
Furthermore, this equation represents a global critical condition for 
the equilibrium problem from the solitonic point of view, as showed 
in (3.35). We remember that the critical condition (3.10) is equal 
to the condition on * 

v*xvf = o v(x, y) g n a n rt 
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It's clear at this point that the Throumoulopoulos et al. and Lapenta 
theoretical positions (see (2.45a-b) and (2.46a-b)) about a specific 
transformations set for the free fields are cannot opposite, because 
both these choices satisfy the critical condition (3.8b), which is iden- 
tically verified for the first choice, while it is in its general form in the 
second case because however it remains verified for an appropriate 
choice of the u and v real functions inside fl a n ftp. 
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